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Abstract. Spectroscopic analyses of 7 SMC B-type supergiants and 1 giant have been undertaken using high resolution optical 
data obtained on the VLT with UVES. FASTWIND, a non-LTE, spherical, line-blanketed model atmosphere code was used to 
derive atmospheric and wind parameters of these stars as well as their absolute abundances. Mass-loss rates, derived from 
H a profiles, are in poor agreement with metallicity dependent theoretical predictions. Indeed the wind-momenta of the SMC 
stars appear to be in good agreement with the wind-momentum luminosity relationship (WLR) of Galactic B-type stars, a 
puzzling result given that line-driven wind theory predicts a metallicity dependence. However the galactic stars were analysed 
using unblanketed model atmospheres which may mask any dependence on metallicity. A mean nitrogen enhancement of a 
factor of 14 is observed in the supergiants whilst only an enrichment of a factor of 4 is present in the giant, AV216. Similar 
excesses in nitrogen are observed in O-type dwarfs and supergiants in the same mass range, suggesting that the additional 
nitrogen is produced while the stars are still on the main-sequence. These nitrogen enrichments can be reproduced by current 
stellar evolution models, which include rotationally induced mixing, only if large initial rotational velocities of 300 kms -1 are 
invoked. Such large rotational velocities appear to be inconsistent with observed v sin i distributions for O-type stars and B-type 
supergiants. Hence it is suggested that the currently available stellar evolution models require more efficient mixing for lower 
rotational velocities. 
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pergiants - stars: mass-loss - stars: abundances - stars: evolu- 
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^ 1. Introduction 



Due to its low metal content and close proximity the Small 
Magellanic Cloud (SMC) is a key laboratory in which to study 
massive stars. Its metallicity (Z), or more accurately the abun- 
dance of its iron peak and alpha-processed elements, is ap- 
proximately one fifth solar, making it an excellent template for 
studies of metal poor starbursts and high redshift star-forming 
galaxies (Leitherer et al. 2001, Pettini et al. 2000). Its proxim- 
ity, at a distance modulus of 18.9 (Harries, Hilditch, & Howarth 
2003), makes high resolution studies feasible with 4-8 m tele- 
scopes and since the extinction is low, E(B-V)^0.07, similar 
studies are also feasible in the important ultraviolet regime. 

This combination of circumstances has led to a number of 
detailed studies of O-type stars in the SMC at optical and UV 
wavelengths (Walborn et al. 1995; Puis et al. 1996; Haser et 
al. 1998). Even more recently we have seen the publication 
of work related to a high resolution UV study of a sample of 
O-stars in the SMC (Walborn et al. 2000; Bouret et al. 2003; 
Hillier et al. 2003). In the present paper, and a companion pa- 
per (Evans et al. 2004; hereafter EVANS04), we extend this 
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work into the domain of the B-type stars in the SMC. There are 
important reasons for this extension, among them the fact that 
the B-type stars are important contributors to the UV spectrum 
of continuous starforming and intermediate age (20-50 Myr) 
stellar populations (deMello, Leitherer, & Heckmann 2000). In 
fact the UV spectra of the stars discussed in the present paper 
have been used to constrain the age of the metal poor super 
star cluster NGC1705-1 (Vazquez et al. 2003). However what 
is more relevant to the present paper is that B-type supergiants 
in the SMC provide important information on mass-loss and 
stellar evolution in a low metallicity environment. 

The surface compositions of B-type supergiants are a key 
test of stellar evolution models since mixing and mass-loss 
processes may alter their surface CNO abundances, especially 
that of nitrogen. The SMC is an ideal laboratory in this re- 
spect since its pristine nitrogen abundance is about l/30th to 
l/50th of solar and enhancements are correspondingly easier to 
detect. Enhanced nitrogen abundances have already been ob- 
served in Galactic and Magellanic Cloud B-type supergiants 
indicating some mixing of CNO-cycled material to the stel- 
lar surface (Gies & Lambert 1992; Lennon et al. 1991, 1996, 
1997, 2003; Fitzpatrick & Bohannan 1993; McErlean, Lennon, 
& Dufton 1999; Dufton et al. 2000). Moreover, Lyubimkov 
(1991) and later Gies & Lambert (1992) found that Galactic 
main-sequence B-type stars also show such N enhancements in 
their spectra. Whether the mechanism which causes this mixing 
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of core processed material to the photosphere is the presence of 
a blue-loop evolutionary phase (Schaller et al. 1992) or due to 
stellar rotation (Denissenkov 1994; Talon 1997) is still uncer- 
tain. At present much effort is being put into the development 
of stellar evolution models which include rotationally induced 
mixing (Heger & Langer 2000; Maeder & Meynet 2000, 2001). 
Such models have had success in predicting the correct depen- 
dence of the blue to red supergiant (B/R) ratio on metallicity, 
previously a problem for models without rotation. This ratio 
is an important factor in stellar and galactic evolution models 
as it represents the relative lifetimes of hot and cool core he- 
lium burning phases. Eggenberger, Meynet, & Maeder (2002) 
showed that evolution models including rotationally induced 
mixing (Maeder & Meynet 2000) seem to reproduce the ob- 
served B/R ratio at Galactic, LMC and SMC metallicities, as 
suggested by Langer & Maeder (1995). Lennon, Dufton, & 
Crowley (2003; hereafter LEN03) investigated surface abun- 
dances of a sample of lower luminosity B-type giants in the 
SMC, and discussed these in the context of evolution with rota- 
tion. Whether these evolutionary models can reproduce the en- 
hanced nitrogen abundances and rotational velocities observed 
in supergiants/giants in a low metallicity regime will be ex- 
plored in this paper. 

Blue supergiant stars, by virtue of their low bolometric cor- 
rection and high intrinsic luminosity, are the most accessible 
stars to study at visible wavelengths. The first chemical com- 
positions of individual stars in the Magellanic Clouds were de- 
termined from B- to F-type supergiants by Przybylski (1968, 
1971), Wolf (1972, 1973) and Osmer (1972, 1973), but with 
mixed results. Osmer (1973) concluded, using simple curve of 
growth arguments, that since the metal lines in a B 1.5 and B2.5 
supergiant were approximately a factor of two lower than their 
Milky Way counterparts then their metallicities must be at least 
a factor of four below solar. This is remarkably close to the 
now commonly accepted value of one-fifth solar and its accu- 
racy is probably due in part to the care with which galactic 
standards were chosen as comparison stars. Subsequently blue 
supergiants have provided important insights into the current 
chemical makeup of many galaxies in the Local Group and be- 
yond (Monteverde et al. 1996, 1997, 1998, 2000; Venn et al. 
1995, 1999, 2000, 2001, 2003a, 2003b; McErlean, Lennon, & 
Dufton 1999; Dufton et al. 2000; Smartt et al. 2001; Bresolin 
et al. 2001, 2002a, 2002b; Trundle et al. 2002; Urbaneja et 
al. 2003) using somewhat more sophisticated modeling tech- 
niques. Lennon et al. ( 1991) used plane-parallel non-LTE meth- 
ods to analyse a B-type supergiant in each Cloud, this being fol- 
lowed by a more comprehensive non-LTE study of 48 B-type 
supergiants in the SMC by Dufton et al. (2000). Nevertheless 
as discussed by McErlean, Lennon, & Dufton (1998) the ef- 
fect of microturbulence and the neglect of the stellar wind in 
this plane-parallel approach are significant sources of uncer- 
tainties in their results. With the recent development of stel- 
lar atmosphere codes which treat in a unified manner both the 
photosphere and wind of massive stars, it is now possible to 
include the effects of the stellar wind (Santolaya-Rey, Puis, & 
Herrero 1997; Herrero, Puis, & Najarro 2002; Repolust, Puis, 
& Herrero 2003 - FASTWIND; Hillier & Miller 1998 - CMF- 
GEN; Pauldrach, Hoffmann, & Lennon 2001 - wm-basic). We 



use FASTWIND in the present paper which in addition to deter- 
mining the chemical composition of our sample, enables us to 
derive estimates for their mass-loss rates and wind momenta. 

Kudritzki et al. (1995) introduced the Wind-Momentum 
Luminosity Relationship (WLR) as a method of determining 
distances out to the Virgo and Fornax Clusters through the 
spectral analysis of hot, massive stars. Line-driven wind the- 
ory implies that photon momentum transfer through metal line 
absorption is driving the stellar wind. As the WLR is based 
on line-driven wind theory, it follows that this relationship 
should be dependent on metallicity. Studies of O-type super- 
giants in the SMC and LMC have confirmed this metallicity de- 
pendence. Additionally, analyses of Galactic early-type super- 
giants have shown some dependence on spectral type (Puis et 
al. 1996; Kudritzki et al. 1999). More recent analyses of Bouret 
et al. (2003) and Repolust, Puis, & Herrero (2003) show a dis- 
tinct difference in the wind momenta of O-type supergiants to 
that of giants and dwarfs in both the Galaxy and Magellanic 
Clouds. In this paper we will compare the WLR of our sam- 
ple with that of stars in the Galaxy and with the predictions of 
theory. 

2. Observations and Data Reduction 

High resolution echelle spectra were obtained for eight SMC 
B-type stars covering the spectral range of B0 - B5, using the 
UVES spectrograph on the VLT. The data were taken during 
twilight on a 3 night run in November 2001, see the obser- 
vational details in Tabled UVES is a two arm crossdispersed 
echelle spectrograph, the red arm of which contains a mosaic 
of two chips. EEV CCD-44 2 k*4 k chips are used in both arms 
along with a MIT-LL CCID-20 CCD chip in the red arm which 
reduces the fringing beyond 7000 A. The useful wavelength 
ranges obtained from the three chips were 3750 - 5000 A, 5900 
- 7700 A & 7750 - 9600 A. These CCD chips have a pixel size 
of 15 /im and a slit width of 1.5 arcsec which was used to give 
a resolving power of approximately 20,000. The spectra were 
later binned to a dispersion resolution of 0.2 A per pixel. 

The 2-dimensional CCD frames were transformed to mul- 
tiple order spectra using the UVES pipeline software. This was 
done in an interactive mode through the ECHELLE/UVES con- 
text of MIDAS. Th-Ar arc lamp exposures were used to define 
echelle order positions and wavelength calibrate the spectra. 
A constant bias level and median background were subtracted 
from the flat-fields which were then extracted in pixel-order 
space. The stellar spectra were reduced using both average 
and optimal extraction methods. In the blue region our anal- 
ysis utilised the optimally extracted spectra, however due to 
the magnitude of the fringing in the red region the averaged ex- 
tracted data proved to be more useful. The final merged spectra 
and individual orders were then processed using an automated 
IDL routine to filter out the cosmic rays. 

Normalisations, corrections for radial velocity shifts and 
equivalent width (EW) measurements were carried out in 
DIPSO, the spectrum analysis procedure of STARLINK. The 
measured equivalent widths for the lines used in the abundance 
analysis are given in Table|2] Although the REDUCE/UVES rou- 
tine normalised and merged the echelle orders, these merged 
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Table 1. Observational details. Identification numbers are from Azzopardi & Vigneau (1982; AV#) and Sanduleak (1968; Sk#). 
Spectral types adopted from Lennon (1997). Absolute Magnitudes (M v ) are calculated using V and (B-V) magnitudes from 
Garmany, Conti & Massey (1987; 1 ) and Massey (2002; 2 ) and (B-V) values from Fitzpatrick & Garmany (1990). The adopted 
distance modulus is 18.9 (Harries, Hilditch, & Howarth 2003). The average S/N ratios for the VLT/UVES data are given, but the 
actual S/N varied across the orders. Heliocentric corrected radial velocities are presented for both the VLT/UVES and NTT/EMMI 
datasets. Note there is no significant variation in radial velocities from the two datasets and thus no obvious indication of binary 
stars amongst the sample. In addition we note that the quoted v sin i values are upper limits and are dominated by macroturbulence 
as discussed in Sect. 2. 



STAR Alternative Spectral V B-V M„ S/N v Jsr vshii 

ID Type B R UVES EMMI (kms -1 ) 



AV215 


Sk76 


BNO la 


12.69 2 


-0.09 


AV104 




B0.5 la 


13.17 2 


-0.16 


AV216 




Bl III 


14.22 2 


-0.13 




Skl91 


B1.5 la 


11.86 1 


-0.04 


AV210 


Sk73 


B1.5 la 


12.60 2 


-0.02 


AV18 


Skl3 


B2 la 


12.46 1 


0.03 


AV362 


Ski 14 


B3 la 


11.36 1 


-0.03 


AV22 


Skl5 


B5 la 


12.25 1 


-0.10 



-6.61 


120 


110 


154 ± 10 


159 ± 15 


91 


-5.82 


150 


120 


163 ± 12 


165 ± 28 


80 


-5.08 


120 


100 


203 ± 14 




77 


-7.41 


170 


140 


130 ± 24 


134 ± 17 


96 


-6.73 


180 


140 


173 ± 8 


192 ± 12 


65 


-7.00 


170 


120 


148 ± 10 


127 ± 13 


49 


-7.82 


190 


100 


208 ± 14 


214 ± 10 


51 


-6.62 


140 


120 


139 ± 12 


109 ± 16 


46 



spectra were checked manually. Low order polynomials were 
fitted to the continuum of each order for normalisation, the 
orders were then merged using appropriate weightings. In 
most cases, the automatically and manually reduced data were 
in good agreement. Nevertheless in measuring the equivalent 
widths of lines in the overlapping regions of different orders 
and in analysing the Balmer lines, more weight was given to 
the manually reduced spectra. 

The metal and diffuse helium absorption lines were fit using 
a non-linear least squares technique with a Gaussian profile. 
For isolated, well observed lines the full-width-half-maximum 
(FWHM) and position were allowed to vary. Blended lines such 
as the C II 4267 A and Mg II 448 1 A were treated as multiple 
Gaussians with fixed wavelength separations and adopting the 
mean FWHMs. For well observed lines error estimates of 10 
% are appropriate, however blended and weak lines have errors 
of 20 % or more. The signal-to-noise ratio (S/N) obtained was 
measured throughout the spectral range of each arm in regions 
of continuum with no observed spectral lines. In the blue region 
S/N is typically 150 or higher and in the red all spectra have 
a S/N greater than 100. Radial velocity shifts were measured 
from the positions of the well observed lines and are presented 
in Table[2along with the mean S/N ratios. 

Additionally, medium resolution spectra of the seven su- 
pergiants were previously obtained using the NTT telescope in 
remote observing mode with the ESO Multi Mode Instrument 
(EMMI). The reduction methods of this dataset were discussed 
in Lennon (1997) and all equivalent widths were remeasured 
here using the technique described above. Each object was 
observed for between 10 - 15 mins. Depending on its visual 
magnitude and the weather conditions at the time, these expo- 
sure times gave signal-to-noise ratios greater than 60 in the red 
and 90 in the blue. The equivalent widths measured from this 
dataset agree with the UVES dataset to within 5 - 20 % depend- 
ing on the quality of the spectra. Comparison of radial velocity 
measurements from the EMMI and UVES datasets show no signs 



of binaries amongst the stellar selection (see Table^J. However 
from STIS/HST target acquisition images, EVANS04 noticed a 
companion to AV216, the faintness of which suggests that it is 
a main-sequence star. Indeed some asymmetry was observed in 
the H Q and He I 4471 A profiles (see Figs.QJ&HJ. The addi- 
tional dataset was mainly used for a check on our UVES data. 
However due to some nebular contamination of the H Q profiles 
of AV362 & AV22, the EMMI dataset was used to derive the 
mass-loss rate of these two stars. There was also some nebular 
contribution observed in the H Q profile of Skl91, and we have 
taken this into account in fitting the data. 

The projected rotational velocity for each star, vsmi, was 
determined from a selection of seven unblended metal lines by 
fitting theoretical profiles to the observed spectra. The theo- 
retical profiles were convolved with a Gaussian profile to ac- 
count for the instrumental broadening. The profile was then 
convolved with a rotational line profile function, this has been 
described in detail by Rucinski (1990) and more recently by 
Gray (1992). Howarth et al. (1997) and Ryans et al. (2002) 
found that the broadening of spectral features in B-type su- 
pergiants was dominated by macroturbulence. Indeed our pro- 
jected rotational velocity estimates listed in TableJQshow a sim- 
ilar behaviour. Ryans et al. (2002) also estimated that the typi- 
cal projected rotational velocity was of the order 10-20 kms -1 . 
Hence the values for v sin i tabulated in Tablefflshould be con- 
sidered as indicative of the width of the spectral features and as 
upper limits on the projected rotational velocities. 

3. Model Atmospheres 

'Unified model atmosphere' codes such as those mentioned in 
the introduction have an advantage over plane-parallel codes 
(viz. TLUSTY - Hubeny & Lanz 1995). This is due to treat- 
ing the atmosphere in a consistent manner by including a 
smooth transition, from a pseudo-hydrostatic photosphere to 
an outer expanding atmosphere - 'the wind'. We have opted 
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Table 2. Estimated equivalent widths (mA) for SMC B-type supergiants. Uncertainties range from 10 % for the best observed 
features to 20 % for the weaker features. 



Line 


AV215 


AVI 04 


AV216 


Skl91 


AV210 


AV18 


AV362 


AV22 


C II 4267 


25 


52: 


49 


48 


83 


81 


53 


55 


C II 6578 










105 


154 


120 


72 


C II 6580 










116 


119 


70 


77 


N II 3995 


68 


49 


40 


137 


173 


169 


95 


81 


N II 4447 








47 


69 


47 


39 


26 


N II 4601 








50 


59 


49 






N II 4607 








29 


61 


46 






N II 4614 










42 


22 






N II 4621 


31 






33 


51 


25 






N II 4630 


73 


48 




95 


158 


103 


53 


47 


O II 4072 


57 


88 


87 


96 


86 


69 






O II 4076 


44 


121 


87 


141 


106 


88 






On 4317 




40 




95 


72 


61 






O II 4320 




36 




87 


87 


63 






O II 4367 


38 


49 




92 


88 


61 






On 4415 




121 


86 


166 


132 


109 






O II 4417 




31 


70 


71 


81 


4066 






O II 4591 


54 


67 


29 


86 


58 


40 






O II 4596 


58 


50 


43 


62 


43 


33 






O II 4638 




40 




141 


81 


169 






O II 4641 




137 




250 


92 


93 






O II 4662 


55 


69 


50 


151 


207 


167 






Mg II 4481 


80: 


65 


46 


75 


114 


93 


182 


217 


Si II 4128 












30: 


73 


94 


Si II 4131 












70: 


79 


96 


Si III 4553 


168 


161 


94 


266 


264 


199 


56 


58 


Si m 4568 


129 


120 


78 


206 


211 


174 


47 


45 


Si in 4575 


64 


54 


40 


108 


113 


95 


21 




Si IV 41 16 


212 


111 


60 


47 


34 









to use FASTWIND (Fast Analysis of STellar atmospheres with 
WINDs), a spherically symmetric non-LTE code first intro- 
duced by Santolaya-Rey, Puis, & Herrero (1997). FASTWIND 
has been developed to include an approximated treatment of 
metal line -blocking and -blanketing (Herrero, Puis, & Najarro 
2002; Repolust, Puis, & Herrero 2003). 

The main advantage of FASTWIND over more sophisticated 
codes such as CMFGEN (Hillier & Miller 1998) is the short 
computational time it takes to calculate models for a significant 
parameter space. Fast computational time allows for quantita- 
tive spectroscopic analyses of large samples of massive stars 
whilst still accounting for the stellar wind. This is achieved 
by parameterising the temperature structure using the non-LTE 
Hopf function described by Santolaya-Rey, Puis, & Herrero 
(1997) which includes the effects of sphericity and mass-loss: 

q'N(rr R ) = ^(tr) (1) 

tr 

where t/r and q/jy are the spherical generalisation of 
the optical depth, tr, and the plane parallel Hopf function, 
1n(tr). qAr(T fl ) and tIr are defined by 

Qn(tr) -9co + (q a - qoD)exp(-jTR) (2) 



and 

driR = XR {r){^fdr (3) 
r 

where q^, q D and 7 are fit parameters to the run of the Hopf 
function, q^, as a function of Rosseland opacity, (xr)- For the 
models discussed in this paper we have adopted Hopf param- 
eters from the parameterisation of non-LTE line-blanketed (in- 
cluding Iron) TLUSTY models calculated at Queen's University 
Belfast (Ryans, R.S.I. ; private communication). The TLUSTY 
grid covered a temperature range of 15 kK to 33 kK and log- 
arithmic surface gravities from 4.5 dex down to the Eddington 
limit. The abundances in this grid were baseline SMC abun- 
dances (viz. [Fe/H]=6.80 dex). The flux conservation reached 
by using the TLUSTY temperature stratification was typically 
of the order of 1 - 2 %, however for stars with high mass-loss 
rates this was higher at <~ 3 %. 

FASTWIND has a similar philosophy to that of the line for- 
mation code DETAIL (Butler & Giddings 1985), i.e the code is 
data driven. The calculated models for this work include atomic 
data for hydrogen, helium, carbon, nitrogen, oxygen, magne- 
sium and silicon. The atomic models and Stark broadening data 
implemented in the line formation of hydrogen and helium are 
described in full in Santolaya-Rey, Puis, & Herrero (1997). In 



C.Trundle et al.: SMC B-type Supergiants. 



5 



Table 3. Derived atmospheric parameters for SMC B-type supergiants. The microturbulence from both the Si III & O II lines 
are presented, however the silicon microturbulence was adopted in our analysis. A microturbulence of 10 kms -1 was adopted for 
AV22, due to the difficulties in constraining £ from the metal lines in its' spectrum. M evo i are estimated from the stellar evolution 
tracks of Maeder & Meynet (2001), see discussion in Sect. 6.2. The errors quoted are the typical random errors as discussed in 
Sect. 4. 



Star 


Spectral 


T„ff 

1 en 




R* 


M 


ML„,„J 
lvl evoL 






fn 
s u 




Type 


(kK) 


(cgs) 


(R©) 


(M ) 


(M ) 


(kms 




AV215 


BNO la 


27.0 ± 1.0 1 


2.90 ±0.10 


30 


26 


39 


5.63 


12 


17 


AV104 


B0.5 la 


27.5 ± 1.0 1 


3.10 ±0.10 


20 


19 


26 


5.31 


11 


20 


AV216 


Bl III 


26.0 ± 1.5 2 


3.60 ± 0.20 


16 


36 


20 


5.00 





15 


Skl91 


B1.5 la 


22.5 ± 1.5 2 


2.55 ± 0.15 


51 


33 


41 


5.77 


13 


20 


AV210 


B1.5 la 


20.5 ± 1.5 2 


2.40 ±0.15 


40 


15 


27 


5.41 


12 


18 


AV18 


B2 la 


19.0 ±2.0 3 


2.30 ± 0.20 


49 


17 


28 


5.44 


9 


16 


AV362 


B3 la 


14.0 ± 1.5 4 


1.70 ±0.15 


96 


17 


30 


5.50 


10 




AV22 


B5 la 


14.5 ± 1.5 4 


1.90 ±0.15 


53 


8 


19 


5.04 


10* 





T c ff determined from : 1 Si Ill/Si IV ionisation balance & He II fits ; 

2 Si Ill/Si IV ionisation balance ; 3 Fits to Si III lines ; 4 Si Il/Si III ionisation balance. 



the case of silicon, nitrogen, and oxygen, the metal ion popu- 
lations and line profiles were calculated using the atomic data 
of Becker & Butler (1988, 1989, 1990). The atomic data used 
for C II are from Eber & Butler (1988), and the Mg II is from 
Mihalas (1972). In most cases these datasets needed some mod- 
ification to ensure convergence for the temperatures and gravi- 
ties considered here, however tests showed that these had neg- 
ligible impact on the resultant line strengths. 

For Si II we initially used the silicon model atom as de- 
scribed by Becker & Butler (1990). This model essentially uses 
the atomic levels adopted by Kamp (1978), but with updated 
atomic data for bound-bound and bound-free transitions (see 
also Lennon et al. 1986). However comparisons of TLUSTY 
and FASTWIND calculations in the low mass-loss rate limit 
(see below) showed severe discrepancies for the Si II equivalent 
widths. This was traced to the simplified atomic model for this 
ion in the FASTWIND calculations. We therefore replaced our 
12-level Si II model ion with a more complete 34-level model 
ion. This has important consequences for the determination of 
the effective temperatures of B-type stars, in that for the mid 
to late B-type stars (B2 - B9) the derived temperatures will be 
cooler using the new Si II model atom. However it is beyond 
the scope of this paper to include a detailed discussion of this 
and so it will be reported in future work. 

To check the validity of our results from FASTWIND, we 
have carried out preliminary tests with the plane-parallel code 
TLUSTY and unified model atmosphere code, CMFGEN (Hillier 
& Miller 1998). Choosing a star with a thin-wind, AV104, and 
a thick-wind, AV215, we analysed these stars with TLUSTY and 
CMFGEN, respectively. No major discrepancies were found in 
these analysis with any differences being within the estimated 
errors. This is reassuring given the different atomic data and 
more importantly the various physical assumptions each code 
adopts (viz. degree of line-blanketing, treatment of wind, cal- 
culation of formal solution etc). 

Finally it should be noted that some of the projected rota- 
tional velocities listed in Tableware relatively large. Therefore, 



an approach allowing for the variation of atmospheric structure 
over the stellar surface (as used for example by Howarth & 
Smith 2001) might be warranted. However as discussed above, 
the actual rotational velocities are likely to be considerably 
smaller and the assumption of spatial homogeneity should be 
appropriate. 



4. Photospheric and Wind Analysis 

4.1. Stellar parameters 

The atmospheric and wind parameters of our stellar sample are 
presented in Tables [3] & IH respectively. The methods for de- 
termining these parameters are described in the following sec- 
tion; however one should note that they are not derived inde- 
pendently but through an iterative process. 

In order to allow an initial assessment of the parameter 
space appropriate for each of the SMC stars, a coarse grid of 
FASTWIND models were calculated to cover the spectral range 
B0 - B5. For each spectral bin, we adopted one effective tem- 
perature and a range of surface gravities, spanning 0.5 dex. 
The microturbulence adopted was 10 kms -1 as appropriate for 
B-type supergiants, if derived from the silicon spectrum (see 
Sect. 4.1.3). The stellar radius was chosen based on the effec- 
tive temperature of the model as well as bolometric corrections 
and colors taken from Fitzpatrick & Garmany (1990). The (3- 
parameter was fixed at 1 and terminal velocities were taken 
from Kudritzki & Puis (2000), where a calibration of terminal 
velocities with spectral type for Galactic stars is presented. Five 
mass-loss rates for each spectral bin were adopted, covering 
thin to thick winds. To obtain initial estimates of the parame- 
ters for each star, the appropriate set of models for its spectral 
type was chosen and a comparison with the observed data gave 
an approximate mass-loss rate and surface gravity. Once initial 
values were defined, finer, smaller grids were run which were 
appropriate to each individual star. 
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Fig. 1. Profile fits to spectral region 4060 - 4150 A. Observed and theoretical data are shown as solid (-) and dotted (• ■ ■) lines, 
respectively. Note the fits to the wings of the Hs line which were used to fix log g and the varying strengths of the Si II & IV lines 
that were the temperature indicators. Also note the absence of the O II multiplet at 4070 A in the coolest stars AV362 & AV22. 
O II is not observed in the spectra of such cool stars, preventing the determination of £o- 



4.1 .1 . Effective temperature, T e g 

In early B-type stars, the silicon lines are the spectral features 
most sensitive to changes in effective temperature, T e ff (see 
Figs. [0&|2] for the variation of the silicon ions line strength 
with spectral type). The standard method of determining T c ff 
in B-type stars is to search for an ionisation balance between 
consecutive stages of silicon ionisation (Kilian et al. 1991). 
Although these lines are sensitive to the silicon abundance, by 
considering the ratio of the EW's for neighbouring ionisation 
stages we eliminate, to first order, this dependence on abun- 
dance. Initially T c g was constrained by fitting theoretical pro- 
files to the Si lines, this was then followed by a comparison 
between EW's of theoretical and observed profiles. The error 
in this method due to observational constraints and random er- 
rors is typically ±1.5 kK, with a few exceptions as discussed 
below. 



For the hotter stars of our sample (i.e. B0 - B 1 .5) the promi- 
nent ionisation stages of silicon are Si III & IV. In the case of 
AV104 & AV215, the He II lines at 4200 and 4542 A were 
used as further constraints on the temperature, assuming a nor- 
mal helium abundance. Given this additional check on T c ff the 
error in these measurements is ± 1.0 kK. In AV362 & AV22, 
the coolest objects in our sample, Si II & III are the dominant 
silicon ions with Si IV absent from the spectra. Unfortunately 
neither the Si II nor Si IV features are very strong in B2-type 
stars at low metallicities such as that of the SMC, in fact we 
only observe noise features in the region of these lines in AVI 8. 
Although we could not accurately fix an ionisation balance, the 
absence of Si II & IV features provided lower and upper limits, 
respectively, for the temperature of AVI 8. Thus, the error we 
adopt in our temperature estimate for AVI 8 is slightly larger 
than for the rest of our sample (± 2.0 kK). 
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Fig. 2. Profile fits to spectral region 4460 - 4600 A. Observed and theoretical data are shown as solid (-) and dotted (• • •) lines, 
respectively. Note the varying strengths of the Si III lines, which were used as temperature indicators. Also note the presence of 
the He II line at 4542 A in AV215 and AVI 04, which was used as an additional temperature diagnostic. The asymmetry in the 
He I 4471 line of AV216 suggests some multiplicity. 



The inclusion of sphericity and line-blanketing in the anal- 
ysis of O-type dwarfs and supergiants has introduced a new 
temperature scale for O-type stars (Martins, Schaerer, & Hillier 
2002; Herrero, Puis, & Najarro 2002; Crowther et al. 2002; 
Repolust, Puis & Herrero 2003). Compared to plane-parallel, 
unblanketed models the temperatures derived for O-type stars 
are up to 8.0 kK cooler, when line-blanketed, spherical mod- 
els are applied. Herrero, Puis, & Najarro (2002) noted that this 
discrepancy decreased to 2.0 kK at the 09.5 la spectral type. 
For the B0 - B2 -type stars in our sample no notable difference 
was observed in our temperature estimates when compared to 
the calibration of Dufton et al. (2000) for B-type supergiants 
(this was also seen in Repolust, Puis & Herrero 2003). However 
for the B3 and B5 star, our temperature estimates were much 
lower than their temperature scale. This is probably a result 
of our improved Si II model, which for a given temperature, 
weakened the Si II line strength. Additionally these two cooler 



stars have morphologically different spectra from which their 
spectral types have been derived, however they have very sim- 
ilar temperatures. The fact that the temperatures of these stars 
are indistinguishable between the two spectral types, is partly 
due to the lower luminosity and hence higher gravity of AV22; 
this temperature 'degeneracy' is also observed between AV215 
and AV104 (B0 & B0.5 stars). An additional factor that may 
be coming into play here is that AV215 and AV362 both have 
very high mass-loss rates. Herrero, Puis, & Najarro (2002) and 
Evans et al. (2003) observed that for extreme supergiants the 
effective temperature derived is much lower than for normal 
supergiants. Voels et al. (1989) also observed a reduction of 
2.0 kK in the temperature of a Cam, an 09.5 la type star, with 
the inclusion of the stellar wind in the analysis. They attributed 
this to 'wind-blanketing', the backscattering of radiation onto 
the photosphere from the stellar wind (see also Repolust, Puis 
& Herrero 2003). 
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4.1 .2. Surface gravity, log g 

The logarithmic surface gravity, log g, is determined in an it- 
erative process with T c ff. The sensitivity of the Balmer line 
wings to gravity provide an excellent tool for deriving log g. 
As H Q and H^ are formed higher up in the atmosphere than H 7 
and H$, these lines are more affected by the wind. Therefore 
the latter lines were used to determine the surface gravity with 
the fitting uncertainty being typically ± 0.05 dex. However the 
uncertainty in the effective temperature estimates induces an 
additional though correlated uncertainty in log g. These error 
estimates are presented in Table [3] along with the logarithmic 
surface gravities, while the theoretical and observed H$ pro- 
files are shown in Fig. Q These measured log g values are in 
fact only the effective gravity, which is reduced due to the cen- 
trifugal acceleration by (v sin i) 2 /R*- The correction to the de- 
rived surface gravities due to the centrifugal force is less than 
0.03 dex for these stars as a result of the low projected rota- 
tional velocities and large stellar radii. As discussed in Sect. 2, 
the measured v sin i value is an upper limit to the projected ro- 
tational velocity and so the correction to the surface gravity for 
these B-type supergiants is negligible. In addition, we note that 
the gravity for AV216 may be overestimated due to the com- 
panion star, which is possibly a main-sequence star (see Sect. 
2) and would therefore have broader Balmer wings. 

4.1 .3. Microturbulence, £ 

Microturbulence, £, is introduced to ensure no variation of 
abundance with line strength from a particular ionic species. 
In early B-type stars there is a sufficient number of O II and Si 
III lines covering a range of line strengths from which to deter- 
mine microturbulence. However, the estimates from different 
metal species can differ by up to 10 kms -1 (see Vrancken et al. 
1998, 2000; McErlean, Lennon & Dufton, 1999; Trundle et al. 
2002). This effect was again observed in our sample where the 
microturbulence deduced from the O II lines, £o, were on aver- 
age 7 kms -1 higher than those derived from the Si III lines, £si- 
Whilst the oxygen lines have a significant range of EW's from 
which to determine the microturbulence, this EW range is not 
present in one particular multiplet but over a number of multi- 
plets. Therefore the scatter in oxygen abundances from one line 
to another is not simply due to microturbulence but to uncer- 
tainties in the modelling (e.g. errors in the oscillator strengths 
or departure coefficients between different multiplets). This is 
not the case for the Si III multiplet which has three lines at close 
wavelength intervals due to the small difference in the energy 
levels involved in these transitions, as such they are effected by 
the same radiation field. Therefore any errors in the oscillator 
strengths or departure coefficients would be similar for each 
line. The microturbulent velocities from both elements are pre- 
sented in Table although we adopt the microturbulence de- 
rived from the Si III triplet at 4560 A (4553, 4568, 4575 A). 
Due to the low T c ff of AV362 and AV22 no O II lines were ob- 
served and so could not be derived (see Fig.Q. In addition, 
for AV22 we had difficulty constraining the silicon microtur- 
bulence as the Si III 4575 A line was not visible in the spec- 
trum and thus there was not a significant enough range in line 



strength to determine £gi. Therefore, we adopted a microturbu- 
lence of 10 kms -1 . A discussion of the effect of microturbu- 
lence on the abundances of our supergiant sample is given in 
Sect. 5. 

In previous non-LTE analyses of supergiants the microtur- 
bulences determined have been significantly lower than val- 
ues derived from LTE analyses, which tended to be above the 
speed of sound (Becker & Butler 1989; Gies & Lambert 1992; 
McErlean et al. 1999; Vrancken et al. 2000; Trundle et al. 
2002). Although for main-sequence B-type stars the value of £ 
is typically lower than 5 kms -1 , in supergiants the microturbu- 
lence is still quite high even with the inclusion of non-LTE ef- 
fects. It was suggested by Lennon et al. (1991) and investigated 
by Kudritzki et al. (1992) that the microturbulence may be re- 
duced for supergiants by the consideration of the wind outflow 
in these stars. However, it is evident from the values in Table|3] 
that microturbulence is still needed in analyses which include 
the non-LTE effects and winds of supergiants. In fact the inclu- 
sion of spherical symmetry and treatment of the stellar wind in 
our model atmospheres doesn't appear to have any effect on the 
magnitude of microturbulence. From a plane-parallel non-LTE 
analysis of a group of galactic B-type supergiants, McErlean, 
Lennon, & Dufton (1999) found an average microturbulence of 
1 1 kms -1 and this is consistent with our values. For AV216, the 
giant in our sample, the Si III lines were compatible with £ = 
kms -1 , while the O II lines implied a £ = 15 kms -1 . We note 
that adopting the latter value would imply a higher effective 
temperature estimate for the Si Ill/iv ionisation equilibrium. In 
turn this would lead to He II line strengths incompatible with 
the observed helium spectrum of AV216. 

4.1 .4. Radius, luminosity, and mass 

Apparent V-magnitudes were taken, where possible, from the 
UBVR CCD survey of Massey (2002). For the brighter tar- 
gets which were not covered in this survey, m^ was taken from 
Garmany, Conti, & Massey (1987). Absolute magnitudes were 
than deduced using (B-V)o from the calibration of Fitzpatrick 
& Garmany (1990) and adopting an SMC distance modulus 
(DM) of 18.9 (Harries, Hilditch, & Howarth 2003). As the in- 
trinsic colours of Fitzpatrick & Garmany (1990) were derived 
for the LMC they are probably not the most suitable values to 
apply to the SMC. Yet given the lack of any calibration for the 
SMC and the errors involved in these calibrations we find they 
are the only appropriate choice. 

The stellar radius was determined using a method described 
by Kudritzki (1980) 

5 log = 29.58 - (M v - V) (4) 
where, 

V = -2.5 log/ F x S° v (\)d\ (5) 
Jo 

F\ is the stellar flux from the final theoretical model of each 
star (units: erg s -1 cm -2 A -1 ) and S°(A) is the V-filter re- 
sponse function for air mass zero from Matthews & Sandage 
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Fig. 3. H a profile fits for the target stars. Shown are the best fits (red line) using the parameters presented in Table|3]&E] The 
green and blue lines represent a change in mass-loss rate of ± 15 %. However, for AV210, the best fit adopting a /3-parameter 
of 3 (red line) and 1 (green line) are plotted. Note the asymmetry in the H Q profile of AV216, this suggests that the companion 
star is contaminating AV216's spectrum. For AV18, Skl91, AV362 & AV22 the best fit was chosen by fitting the emission peak. 
The latter three have nebular contributions in the H Q profile which were avoided in the fitting procedure. Note how for AV362 & 
AV22 the models have a large absorption trough which is not present in the observed spectrum. 



(1963). The stellar luminosities and masses were then calcu- 
lated using the effective temperatures, gravities, and radii de- 
termined from our analysis. 

The error on the stellar radius is dominated by the un- 
certainty in our distance estimate. The distance to the SMC 
has been determined from many different astrophysical objects 
(viz. Cepheids, eclipsing binaries & RC stars) and has been 
reviewed thoroughly in Harries, Hilditch, & Howarth (2003). 
From the various estimates the SMC lies at a DM between 18.7- 
19.1, therefore we adopt an error of ± 0.2. With an additional 
error of 0. 1 from the adopted apparent magnitudes, this infers 
an error of ± 0.22 on the absolute magnitude estimates giv- 
ing a 10 % uncertainty in the stellar radii. This uncertainty in 
the radius introduces larger errors in the luminosity and mass 
estimates, on the order of 15 % and 30 %, respectively. A sys- 



tematic error is also introduced to the derived mass-loss rates 
due to the uncertainty in the adopted distance to the SMC, this 
is discussed further in Sect. 4.2.2. 



4.2. Wind parameters 

In addition to the stellar parameters which describe their pho- 
tospheric conditions, the spectra provide information on their 
stellar winds. The stellar wind is parameterised by the mass- 
loss rate (M), the exponent of the velocity law (/3-parameter), 
and the terminal velocity (Vqo). The former two parameters are 
derived here from the H Q profile and the latter from UV P- 
Cygni profiles. 
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Table 4. Derived wind parameters for SMC B-type supergiants. Escape velocities are taken from EVANS04 for AV215, Skl91, 
AV210 and AVI 8. For the other four they are calculated in this work using the procedure outlined in EVANS04. Terminal 
velocities, (v^), for AV215, Skl91, AV210, and AV18 are taken from EVANS04, as derived from the SEI method. For the rest 
of the sample, v^ is calculated by adopting a Voo/v csc ratio of 2.65 for stars with T e g > 21.0 kK and 1.40 for T c g < 21.0 kK. 
The errors quoted are the typical random errors as discussed in Sect. 4. 



Star 


Spectral 


M 


Voo 


Vesc 


P 


log(D M OM) 




Type 


(10" 6 M yr" 1 ) 


(kms" 1 ) 


(kms -1 ) 




(cgs) 


AV215 


BNO la 


1.35 ±0.14 


1400 ± 100 1 


447 1 


1.4 ±0.2 


28.81 ±0.10 


AV104 


B0.5 la 


0.40 ± 0.05 


1340 ± 270 2 


506 2 


1.0 ±0.3 


28.18 ±0.15 


AV216 


Bl III 


0.12 ±0.02 


2410 ± 480 2 


908 2 


1.4 ±0.3 


27.86 ±0.15 


Skl91 


B1.5 la 


0.68 ± 0.07 


425 ± 50 1 


374 1 


2.0 ± 0.2 


28.11 ±0.10 


AV210 


B1.5 la 


0.20 ± 0.05 


750 ± 100 1 


290 1 


3.0 ±0.3 


27.52 ±0.15 


AV18 


B2 la 


0.23 ± 0.03 


325 ± 50 1 


284 1 


3.0 ±0.2 


27.66 ±0.10 


AV362 


B3 la 


0.80 ± 0.08 


270 ±50 2 


192 2 


1.0 ±0.3 


28.12 ±0.15 


AV22 


B5 la 


0.23 ± 0.02 


280 ±60 2 


198 2 


1.0 ±0.3 


27.47 ±0.15 



Terminal and escape velocities as measured in 1 EVANS04 and estimated in 2 this work. 



4.2.1. Terminal velocity, 

Terminal velocities of AV215, Skl91, AV210 & AV18 have 
been presented in EVANS04, where a more detailed discussion 
can be found. Briefly, v^ was derived from the UV P Cygni 
lines using the method described by Haser et al. (1995). This 
procedure utilises the Sobolev Exact Integration (SEI) method 
introduced by Lamers, Cerruti-Sola, & Perinotto (1987; see 
also Hamann 1981). For early B-type stars the doublets Si 
IV 1394, 1403 and C IV 1548, 1551 display the most promi- 
nent wind features and thus were the diagnostic lines used by 
EVANS04, see the fits to these profiles in Fig. 2 of EVANS04. 

Unfortunately four stars in our sample showed no signif- 
icant wind features in the aforementioned UV lines; AV104, 
AV216, AV362 & AV22 (EVANS04). Therefore terminal ve- 
locities for these stars were estimated by an alternative method. 
Radiation driven wind theory (Castor, Abbott, & Klein 1975; 
Abbott 1978) predicts that terminal and escape velocities scale 
together as a function of a, such that foo/^osc — d),if 

the finite cone-angle subtended by the stellar disk is accounted 
for. a is dependent on the effective temperature and is de- 
fined as the effective value of the force multiplier parameter, a, 
defining the exponent of the line strength distribution function. 
Lamers, Snow, & Lindholm (1995) found, for stars with tem- 
peratures > 21 kK, that the Vqo/Vosc ratio was approximately 
2.6 but that at lower temperatures (< 20 kK) this dropped to 
1.3. Kudritzki & Puis (2000) found, from a comparison of the 
Lamers sample of 68 stars with an additional set of stars (~ 
200) from Prinja & Massa (1998), that above 21 kK this ratio 
is 2.65 and below 21 kK it is 1.4. 

In EVANS04 a comparison of Voq/Vcsc derived from the 
Lamers and Prinja & Massa samples discussed above with 
those derived from SMC OB-type supergiants, showed no con- 
vincing variation in this relationship due to the different metal- 
licity environments. The average Voq/Vcsc from the SMC OB- 
type supergiants with T c ff > 21 kK is within the scatter of the 
Galactic star results. Therefore in order to derive the terminal 
velocities of AVI 04, AV216, AV362 & AV22 we have adopted 



the Voo/Vesc ratio of 2.65 for the former two stars and 1.4 for 
the two latter, cooler stars. The escape velocities were derived 
using the same procedure as described in EVANS04 (see equa- 
tions 1 & 2), where v esc was calculated from the ratio of stellar 
mass to radius and corrected for the radiative acceleration due 
to electron scattering. We should note that in all three sam- 
ples (Lamers, Snow, & Lindholm 1995; Prinja & Massa 1998; 
EVANS04) there is a large scatter of approximately 20 % in the 
Voo/Vesc at a given temperature. This is reflected in the the error 
estimates for Vqo given in Table |4]for the four stars where we 
have used this method. 

4.2.2. Mass-loss rate & /3-parameter 

The mass-loss rate is derived mainly from the morphology of 
the H a profile, with less weight given to Hp profile. This pro- 
cedure has been described in detail in Puis et al. (1996) and 
Kudritzki et al. (1999). The /3-parameter effects the central 
emission core of the H Q profile and also its FWHM. Therefore 
for thick winds (i.e when there is emission in the H a profile), 
the /3-parameter and mass-loss rate can be obtained simulta- 
neously. The derived mass-loss rates are presented in Table |4] 
the uncertainties quoted are the random errors involved in the 
fitting process which are typically in the range 10 - 15 %. In ad- 
dition we should consider the systematic error incurred through 
the adopted distance modulus. The uncertainty in DM follows 
through to the mass-loss rate from the 10 % uncertainty it in- 
duces on the stellar radius. If we assume that we get identical 
profiles as long as Q is constant (i.e. that we are fitting Q and 

3/2 

not M), where Q = MJRJ , then the systematic error in the de- 
rived mass-loss rates is ± 15 % (Puis et al. 1996). Fits to the 
H Q profiles are presented in Fig. [3] where we show the best 
fits (red line) along with models representing a change in M 
of ± 15 % (green and blue lines). Some nebular contamination 
was detected in Skl91, AV362 & AV22 from comparisons be- 
tween the UVES and EMMI datasets. After careful extraction of 
the nebular contribution to the H Q profiles the EMMI data were 
considered more reliable for AV362 & AV22 whilst for Skl91 
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Table 5. Derived non-LTE absolute abundances for SMC B-type supergiants. Abundances are given as log [N(X)/N(H)] + 12, 
where X represents the appropriate element. Numbers in the parenthesis represent the number of lines used in the abundance 
analysis. The errors represent the standard deviation of the mean and account for systematic errors as discussed in Sect. 5. 



Star 


C II n 


N II n 


On n 


Mg II n 


Si II 


n 


Si III n 


Si IV n 


AV215 


6.91 ±0.10(1) 


7.96 ±0.12 (3) 


7.97 ± 0.14 (6) 


<7.20 (1) 






7.10 ±0.15 (3) 


7.14 ±0.12(1) 


AVI 04 


<7.00 (1) 


7.41 ±0.15 (2) 


8.05 ± 0.32 (12) 


6.90 ± 0.05 (1) 






6.70 ±0.11 (3) 


6.67 ± 0.22 (1) 


AV216 


6.76 ±0.16(1) 


7.14 ±0.20(1) 


8.22 ±0.39 (7) 


6.64 ± 0.08 (1) 






6.66 ± 0.22 (3) 


<6.71 (1) 


Skl91 


6.89 ±0.10(1) 


7.63 ±0.19 (6) 


8.20 ±0.18 (12) 


6.98 ±0.06(1) 






6.75 ±0.16(3) 


6.50 ± 0.25 (1) 


AV210 


6.93 ±0.14 (3) 


7.60 ±0.12 (7) 


8.18 ±0.18 (12) 


6.89 ± 0.07 (1) 






6.69 ±0.14 (3) 


6.70 ±0.11 (1) 


AV18 


6.98 ±0.18 (3) 


7.50 ± 0.30 (7) 


8.26 ± 0.21 (12) 


6.75 ±0.10(1) 


<7.03 


(2) 


6.81 ±0.11 (3) 




AV362 


7.12 ±0.48 (3) 


8.22 ± 0.25 (3) 




6.72 ± 0.06 (1) 


6.60 ±0.13 (2) 


6.66 ±0.19 (3) 




AV22 


7.07 ± 0.27 (3) 


7.92 ±0.16 (3) 




6.83 ±0.09(1) 


6.75 ±0.16(2) 


6.61 ±0.15(2) 





the UVES data appeared less contaminated. In fitting the H Q 
line of these stars any obvious nebular component was avoided 
(see Fig.|3}. 

For those stars with thin winds, (AV104 & AV216), Ha was 
used as a constraint on the /3-parameter, whilst the core depth 
of the H Q line constrained M (see Fig. [3}. For too high a (3- 
parameter, an emission peak appears in the H^ theoretical pro- 
file which is not observed. This method provides an upper limit 
for the /3-parameter. 

In the case of AV215 (BNO la), the profile shape of the 
H Q line provided a means of determining both M and (3 simul- 
taneously. However, the models could not simultaneously re- 
produce the depth and width of the absorption core. Kudritzki 
et al. (1999) noticed this problem in a Bl lab-type supergiant, 
HD 1 3854, suggesting it may be a lack of line-blanketing which 
was affecting the stratification of the source function of H Q 
throughout the wind. This is not supported by our analysis as it 
has been carried out with a version of FASTWIND that includes 
an approximation of the line-blanketing and blocking. This in- 
ability to reproduce the depth of the blueward absorption core 
is also seen for AVI 8 (B2 la). This star is similar to the B2 la- 
type stars HD 41 1 17 & HD 14143 discussed by Kudritzki et al. 
(1999), in that its emission peak could be fitted giving M and 
j3, whilst the absorption core could not be matched. 

The H a profiles of the two B1.5 la type stars illustrate the 
extreme differences between the winds of stars with the same 
spectral and luminosity classification; whilst Skl91 is in emis- 
sion, AV210 is in absorption. In Fig. [3] we show two fits to 
AV210, one with /? = 1 and the other with /? = 3. Although, 
these are extremely different profiles the mass-loss rate only 
varies by 25 %. We have adopted the parameters from the j3 = 
3 case as the observed Hp line is reproduced and the H a ab- 
sorption core depth, though not the profile shape, can also be 
matched. 

In the coolest stars, AV362 & AV22, M and (3 were derived 
from the height and width of the H Q emission. Note in Fig. [3] 
that the code predicts large absorption troughs not observed 
in the actual data, yet they are an inevitable consequence of 
the Lyman continuum and the Lyman resonance line becom- 
ing optically thick at lower temperatures and large mass-loss 
rates (e.g. Kudritzki & Puis 2000). Although there is probably 



some nebular contamination in the H Q profile of both stars, it 
is unlikely to have masked such a large absorption core. One 
can speculate that this problem might be resolved by includ- 
ing the effects of clumping in the analysis, since the mass-loss 
rates derived here would be overestimated and the Lyman opti- 
cal depths might be reduced. We will come back to this again in 
Sect. 6.3 however a detailed discussion of the effects of clump- 
ing is beyond the scope of this work. 



5. Chemical Abundance Analysis 

In previous investigations of the atmospheres of super- 
giants, their has been clear evidence of chemical processing 
(McErlean, Lennon, & Dufton 1999). The photospheres of 
these luminous stars have chemical patterns which indicate 
contamination of the stellar surface by the products of CN- 
cycle burning and thus we have carried out an abundance anal- 
ysis on this stellar sample to study this further. Adopting the 
atmospheric parameters presented in the previous section we 
derived the chemical composition of our stellar sample. The 
results of the absolute abundance analysis are presented in 
Table|5] We have also compared our mean abundances to those 
derived for B-type giants and main-sequence stars in NGC330, 
an unprocessed B-type star AV304 and SMC H II regions, see 
Tabled 

The most dominant uncertainty on B-type supergiant abun- 
dances is the adopted micro turbulence. As discussed in Sect. 
4.1.3 each element implies a different value for the microtur- 
bulence (see Table[3}. As we have adopted the microturbulence 
from the Si III lines, £gj, it is necessary to explore the effect 
of this choice of microturbulence on the other elements. Firstly 
the uncertainty in £si is ± 2 kms -1 , which introduces an error 
of ~ 0.05 for weak lines and typically ± 0.10 for strong lines. 
For example the carbon abundance, which depends mainly on 
the weak C II feature at 4267 A only varies by ± 0.02 for the 
B0 - B1.5 stars. As this line becomes stronger for the B2 - 
B5 -type stars the effect of microturbulence is larger increas- 
ing to ± 0.07. The largest effects of the uncertainty in £si are 
on the abundances of Si IV in AV215 & AVI 04 and nitrogen 
in AV362 & AV22. If we adopt the microturbulence from the 
oxygen lines, there is a decrease in the abundance of carbon, ni- 
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Table 6. Comparison of mean abundances for the SMC B-type supergiants with previous studies of abundances in the SMC. 
Along with our mean B-type supergiant abundances as derived from Table|3]we show the mean abundances of A-type supergiants 
from Venn et al. (1999, 2003b; Si corrected for non-LTE effects.), SMC main-sequence B-type star AV304 (Rolleston et al. 2003; 
corrected for non-LTE effects), NGC330 cluster main-sequence B-type stars (LEN03; corrected for non-LTE effects), and SMC 
H II regions from Kurt et al. (1999). The errors represent the standard deviation of the mean. 





Supergiants 


B stars 


Hll 


Metal 


This work 


A-type 


NGC330 


AV304 


regions 


C 


7.30 ± 0.09 




7.26 ±0.15 


7.41 ±0.18 


7.53 ± 0.06 


N 


7.67 ± 0.27 


7.52 ± 0.10 


7.51 ±0.18 


6.55 ±0.01 


6.59 ± 0.08 


O 


8.15 ± 0.07 


8.14 ±0.06 


7.98 ±0.13 


8.16 ±0.33 


8.05 ± 0.05 


Mg 


6.78 ± 0.16 


6.83 ± 0.08 


6.59 ±0.14 


6.73 




Si 


6.74 ±0.11 


6.92 ± 0.15 


6.58 ± 0.32 


6.74 ± 0.03 


6.70 ± 0.20 



trogen, and magnesium of normally less than 0. 1 dex. However 
the silicon and oxygen abundances are decreased by up to 0.2 
dex depending on the degree of change in the microturbulence 
and the strength of the spectral features. 

5. 1. Helium 

Achieving a consistent fit to both the singlet and triplet neutral 
helium lines in a stellar spectrum is a known problem. Voels 
et al. (1989) concluded that the extension of the stellar atmo- 
sphere was causing this effect and could be resolved by the in- 
clusion of sphericity in the model atmospheres. Studies of OB- 
type supergiants which omit microturbulence tend to find very 
high helium fractions, of up to 0.4, in order to get a fit to all He 
I lines (Lennon et al. 1991; Herrero et al. 1992). By including a 
non-zero microturbulence McErlean, Lennon, & Dufton (1998) 
found that a more consistent, although not perfect, fit to the he- 
lium lines could be achieved with a solar or near solar helium 
fraction of 0.09. Due to the difficulty in constraining the he- 
lium abundance from the neutral helium lines we have adopted 
a solar helium fraction of 0.09. No significant variation from 
this helium abundance was noted in our sample, although con- 
sistent fits to the helium spectra were not always achieved. We 
should note that the helium abundance has an indirect effect 
on the temperature and gravity and hence the mass-loss rate. 
McErlean, Lennon, & Dufton (1998) showed that by changing 
the helium fraction from 0.1 to 0.2, the T e ff and log g were 
decreased by 0.5 kK and 0.05 dex, respectively, which would 
change M by less than 10 %. 

5.2. Carbon, nitrogen & oxygen 

The carbon abundances are based primarily on the C II 4267 A 
line, with the 6580 A multiplet being visible for the B2 - B5 
stars (see Fig. [5). From Table [5] we see that the mean carbon 
abundance, from those stars for which it is derived solely from 
the 4267 A line (B0 - B1.5), is ± 0.1 dex lower than that de- 
termined by including the 6580 A multiplet (B2 - B5). In the 
coolest stars the 6580 A multiplet actually gives abundances 
of ~ 0.25 dex higher than those determined from the 4267 A 
line. The 4267 A line is known to be particularly sensitive to 



non-LTE effects (Eber & Butler 1988), resulting in lower abun- 
dances than derived from other carbon lines (viz. 3919, 3921, 
4650 A). This has been found by LEN03 in the non-LTE anal- 
yses of B-type stars in NGC330 and that of an unprocessed 
B-type star AV304 (LEN03; Rolleston et al. 2003). However 
the carbon abundances for the stars in our sample show little 
scatter and have a mean abundance for the sample of 6.96 dex. 

The mean oxygen abundance is in agreement with these 
previous studies, the large errors in Table [5] represent the un- 
certainty caused by adopting £si- No oxygen abundance could 
be determined for the AV362 & AV22, as these are too cool for 
the presence of visible oxygen lines in their spectra. 

The nitrogen abundances for our sample exhibit a large 
scatter with the mean being 7.67 dex, a factor of 13 more over- 
abundant than AV304. 

5.3. a-process elements; magnesium & silicon 

Magnesium abundances are based solely on the Mg II 448 1 A 
doublet but are reasonably consistent throughout our sample. 
For AV215, the magnesium line was weak and blended with a 
He I line, thus providing us with only an upper limit. 

The silicon abundance is very sensitive to the adopted tem- 
perature and this is clear from the range of abundances pre- 
sented in Table [5] However, the mean abundance agrees well 
with those of AV304 and SMC H II regions. 

6. Discussion 

6. 1. Chemical composition of SMC 

Before analysing our samples composition, it is useful to dis- 
cuss the baseline composition of the SMC to which we can 
relate our absolute abundances. Rolleston et al. (2003) showed 
that the main-sequence star AV304 has a chemical composi- 
tion similar to that derived from SMC H II regions. Since this 
star shows no signs of processing in its photosphere it should 
represent the present-day composition of the SMC. However 
this was an LTE analysis and in order to compare our non-LTE 
abundances with those of AV304 we need to take into account 
the non-LTE effects. LEN03 carried out such corrections on 
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AV304 and their NGC330 B-type star sample. In addition they 
applied corrections to the carbon abundance to account for the 
lower abundances derived from the 4267 A line. As mentioned 
in Sect. 5.2 the absolute carbon abundance is effected by the 
problematic line 4267 A, which gives abundances lower than 
those derived by weaker C II lines. Thus a direct comparison of 
the carbon abundance derived from 4267 A and that from H II 
regions is inappropriate. Gies & Lambert (1992) and Vrancken 
et al. (2000), derived carbon abundances from weak lines and 
the 4267 A line, respectively. From a comparison of these stud- 
ies LEN03 found the weak lines gave abundances which were 
0.34 dex higher than those determined from the 4267 A line and 
subsequently applied this correction to the carbon abundances 
of AV304 and the NGC330 stars. We have also applied this cor- 
rection to our derived carbon abundances which rely mainly on 
the 4267 A line. In Table [6] we present the corrected AV304 
abundances that can be used as a baseline SMC chemical com- 
position along with the SMC H II region abundances (Rolleston 
et al. 2003; Kurt et al. 1999). Note the excellent agreement be- 
tween AV304 and the H II regions. 

Also included in Tableware the mean abundances from our 
B-type supergiant sample, from SMC A-type supergiants(Venn 
et al. 1999, 2003b) and the corrected abundances of NGC330 
B-type stars (LEN03). The Si abundance presented in Venn 
et al. (1999) was derived from LTE calculations, corrections 
similar to those carried out for AV304 and the NGC330 stars 
have been applied to this value to account for non-LTE effects. 
Comparing the abundances of the a-process elements (Mg & 
Si) in each set of objects we see, as expected, that they are 
in reasonable agreement indicating no processing of these ele- 
ments. In fact assuming that the galactic standards of LEN03 
and Rolleston et al. (2003) represent the galactic B-type star 
abundances, we find that the a-processed elements in our sam- 
ple are 0.42 dex lower than the galactic objects. This is in agree- 
ment with the 0.44 dex difference found between HR2387 and 
AV304 (Rolleston et al. 2003). 

The most significant result is the variation in the nitrogen 
abundances. Whilst the AV304 nitrogen abundance and that of 
the H II regions agree, those for the B & A-type supergiants 
and the NGC330 stars are a factor of 9 to 14 greater. It is inter- 
esting to note that within the B-type supergiant sample itself, 
a sizeable variation in nitrogen abundance is also present (see 
Table [5}. These enhancements are much greater than any er- 
ror incurred in our analysis and are therefore considered to be 
real. The most likely explanation for such an enhancement of 
nitrogen is that we are seeing the CNO processed material con- 
taminating the chemical composition of the photosphere. If this 
is the case, with such high nitrogen abundances we might ex- 
pect a corresponding but smaller enhancement in helium with 
depletions of carbon and possibly oxygen, keeping the sum of 
CNO nuclei constant. However a comparison of the oxygen 
abundance with that of AV304 shows excellent agreement, with 
no indication of an underabundance present. On the other hand 
carbon appears to be underabundant by 0.11 and 0.23 in com- 
parison to AV304 and the H II regions, respectively. This is 
comparable to the carbon depletions produced in the Maeder 
& Meynet (2001) stellar evolution models from similar nitro- 
gen enhancements (< 0.25 dex). Nevertheless given the un- 
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Fig. 4. HR-diagram, luminosity as a function of temperature. 
Shown are the stellar evolution tracks of Maeder & Meynet 
(2001) for an assumed initial rotational velocity of ( — ; solid 
lines) and 300 kms -1 (- -; dashed lines) at stellar masses of 20, 
25, 40, & 60 M Q . Included are the results from our early and 
mid B-type supergiants (•; *),with errors in luminosity repre- 
senting ± 15 % , SMC O-type supergiants of Hillier et al. (A 
; 2003) and Crowther et al. (o ; 2002), and the SMC O-type 
dwarfs of Bouret et al. (o ; 2003) 



certainties in the carbon abundance and the magnitude of the 
correction applied to carbon due to the 4267 A line, this slight 
underabundance may not be significant. The CNO pattern of 
our B-type supergiant sample is in good agreement with that of 
the A-type supergiants and NGC330 B-type stars i.e. N is en- 
hanced, whilst C & O appear to be of normal SMC metallicity. 

6.2. Evolutionary status 

To help understand the evolutionary status of these B-type su- 
pergiants/giant we have compared our results to the stellar 
evolution tracks of Maeder & Meynet (2001) which include 
rotationally induced mixing (see Figs. |4] |6] & EJ. As noted 
by LEN03 the initial abundances of these models are scaled 
to solar abundances and so by adopting 0.2 Z Q they greatly 
overestimate the initial nitrogen abundance i.e. instead of an 
SMC nitrogen abundance of 6.55, they use one-fifth solar or 
7.3 dex. We have recalibrated these models to SMC nitrogen 
abundances, following the discussion of LEN03, who assumed 
the initial nitrogen abundance has no effect on the amount of 
excess nitrogen produced and simply rescaled it to the SMC 
metallicity. 

The position of our stellar sample on the HR-diagram is 
shown in Fig.|4] overlayed on the rotating (- -) and non-rotating 
( — ) tracks of Maeder & Meynet (2001). Note that all our stars 
lie redward of the main-sequence phase, with the two coolest 
stars having positions consistent with the core-helium burn- 
ing phase (this phase is represented by the bump at log T c g = 
4.2 - 4.3). By interpolating between the rotating stellar evolu- 
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Fig. 5. Comparison of evolutionary and spectroscopic masses 
as a function of luminosity. Included are the results from our 
early and mid B-type supergiants (•; ★) and NGC330 B-type 
stars from Lennon et al. (* ; 2003). The error bars represent the 
error in the surface gravity. Also shown is the uncertainty in 
mass due to the adopted distance modulus (- -). 



tion tracks we estimated the evolutionary masses (M evo i) of the 
sample, which are listed in Tables[3] Figure[5]illustrates the dif- 
ference between M evo i and the derived spectroscopic masses, 
M spec , which were determined from the stellar radii and sur- 
face gravities (see Table |3j. This plot implies a mass discrep- 
ancy where M evo i > M spec , for most of the sample. LEN03 
suggested that this mass discrepancy, could be accounted for 
by the uncertainties in the surface gravity which were consid- 
erably larger in their work (due to low S/N). Moreover the 
spectroscopic masses are dependent on the distance modu- 
lus (DM) which introduces an additional systematic error. In 
Fig.[5]the solid and dashed line represents the zero point of log 
(M evo i/M spec ) and the errors on this value due to the adopted 
DM. The errors in the gravity and distance modulus can ac- 
count for some of this mass discrepancy but there remains an 
inconsistency of 0. 1 - 0.2 dex. We note that one of our stars, 
AV216, has M evo i < M spec , similar to one star from LEN03. 
As discussed in Sect. 4. 1 .2 we may have overestimated the sur- 
face gravity of this star due to the presence of a companion and 
this would, in turn, lead to an overestimation of M spec . 

In Fig. [6] we plot the recalibrated stellar evolution tracks, 
showing the development of the surface nitrogen abundances as 
a function of temperature. We have included the photospheric 
nitrogen abundances of all the SMC OBA-type supergiants and 
dwarfs available in the literature (Hillier et al. 2003; Crowther 
et al. 2002; Bouret et al. 2003; LEN03; Venn et al. 2003b) as 
well as the results from our B-type supergiants. The initial rota- 
tional velocity of the stellar models which can reproduce nitro- 
gen abundances similar to those observed needs to be at least 
300 kms -1 and these models show no significant slowdown by 
the end of the hydrogen-burning stage. Although these rapidly 
rotating models imply significant nitrogen enhancements some 
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Fig. 6. Surface nitrogen abundances as a function of tem- 
perature. Shown are the stellar evolution tracks of Maeder & 
Meynet (2001) for an assumed initial rotational velocity of 300 
kms -1 and various stellar masses. Included are the results from 
our early and mid B-type supergiants (•; *), SMC O-type su- 
pergiants of Hillier et al. (A ; 2003) and Crowther et al. (o ; 
2002), SMC O-type dwarfs of Bouret et al. (o ; 2003), NGC330 
B-type stars of Lennon et al. (* ; 2003), and A-type supergiants 
from Venn et al. (+; 2003b). The error bars illustrate the random 
and systematic errors on the nitrogen abundance. 



of the OBA-type supergiants have abundances up to a factor of 
five larger than predicted. 

Fig.0shows the observed projected rotational velocities of 
the OBA-type supergiants and dwarfs against the predicted ro- 
tational velocities. Note that we are probably overestimating 
the values for the supergiants due to spectral line broadening 
being dominated by macroturbulence (see Sect. 2; Howarth et 
al. 1997; Ryans et al. 2002). The observed velocities are sig- 
nificantly lower than expected from the stellar evolution mod- 
els. In order for these stars to have rotational velocities similar 
to those predicted, the inclination of the stellar rotation axis 
would have to be < 20°. Assuming that the distribution of the 
rotational axis, i, is random, it is highly unlikely that our sam- 
ple of stars have such large rotational velocities and that we are 
observing them all pole on. 

The O-type dwarfs of Bouret et al. (2003) have masses sim- 
ilar to the masses of our sample, as such we will consider them 
to be representative of the progenitors for the B-type super- 
giants. Since these O-type stars show a similar range of nitro- 
gen enhancements to the B-type supergiants, it is reasonable 
to assume that the additional nitrogen is produced while the 
stars are still on the main-sequence. If the nitrogen enhance- 
ment is truly a product of rotationally induced mixing, this 
process needs to be more efficient at lower rotational veloci- 
ties than predicted by Maeder & Meynet (2001), see also the 
discussion of Herrero & Lennon (2002). 

A factor that needs to be carefully considered in stellar evo- 
lution models is the effect of the mass-loss rate. This is incorpo- 
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Fig. 7. Rotational velocities as a function of temperature. 
Shown are the stellar evolution tracks of Maeder & Meynet 
(2001) for an assumed initial rotational velocity of 300 kms -1 
and various stellar masses. Included are the v sini values from 
our early and mid B-type supergiants (•; *), SMC O-type su- 
pergiants of Hillier et al. (A ; 2003) and Crowther et al. (o ; 
2002), SMC O-type dwarfs of Bouret et al. (o ; 2003), NGC330 
B-type stars Lennon et al. (* ; 2003), and A-type supergiants 
from Venn et al. (+; 1999). 



rated in the Maeder & Meynet (2001) models and is responsi- 
ble for the loss of angular-momentum through the stellar wind. 
These models use mass-loss rates from Kudritzki & Puis (2000) 
and are modified by rotation using an expression dependent on 
a, a force multiplier parameter. Maeder & Meynet assume that 
a = 0.6 but as discussed in Sect. 6.3 this quantity is not reli- 
ably established for B-type stars. Additionally Kudritzki & Puis 
(2000) used only unblanketed analyses of the O and B type su- 
pergiants, whilst more recent analyses will give more accurate 
mass-loss rates due to their inclusion of line-blanketing in the 
model atmosphere codes (see Sect. 6.3 for a further discussion). 
Clearly some allowance for the variation of mass-loss rates as 
a function of metallicity should also be included in the stellar 
evolution models. 

Binary star evolution models by Wellstein, Langer, & 
Braun (2001), can also produce nitrogen enriched blue super- 
giants with similar temperatures and luminosities to the SMC 
supergiants in our sample. Through mass transfer from the pri- 
mary star, the low mass secondary in a binary system becomes 
more massive. These stars may appear to have lower spec- 
troscopic masses than expected from their luminosity. Thus 
this could account for some of the mass discrepancy in our 
stellar sample. One advantage of the binary evolution mod- 
els, is that they can place these blue supergiants on the post 
main-sequence gap which presently can not be achieved with 
single-star models due to their rapid evolution through this part 
of the HR-diagram. Since the primary of the binary system 
will evolve into a low mass helium star by the time the sec- 
ondary is a blue supergiant, it is unlikely that the binary systems 



Fig. 8. Comparison of theoretical and observed mass-loss 
rates. Theoretical predictions are calculated using the metal- 
licity dependent mass-loss recipes of Vink et al. (2001) for Z 
= 0.2 Zq. The solid circles (•) represent the SMC early B-type 
supergiants, whilst the solid stars (*) represent the mid B-type 
supergiants. 



could be detected through radial velocity variations. Wellstein, 
Langer, & Braun (2001) predict radial velocity variations of ~ 
lSkms- 1 . 

6.3. Winds of B-type supergiants 

According to radiative line-driven wind theory, the wind- 
momentum and hence mass-loss rate should be dependent on 
metallicity and may be expressed as a function of the force 
multiplier parameter, a. This parameter is the ratio of the opti- 
cally thick line acceleration to the total line acceleration and the 
metallicity dependence of the mass-loss rate is found to scale 
as 

M~(Z*/Z )^ (6) 

where, a' = a - 5 and 6 is a force multiplier parameter that 
accounts for the ionisation balance in the wind (Kudritzki et 
al. 1989; Puis, Springmann, & Lennon 2000). Theoretically, 
this metallicity dependence is predicted to be in the range 
Z 5 — Z° - 85 for OB-type supergiants (Kudritzki, Pauldrach, 
& Puis 1987; Leitherer, Robert, & Drissen 1992; Vink et al. 
2001). Observationally, a significant enough sample of stars 
has not yet been analysed to get a reliable estimate of a. 
However for three SMC O-type supergiants Puis et al. (1996) 
derived an exponent of 0.52, whilst for 6 O-type dwarfs anal- 
ysed by Bouret et al. (2003) a value of 0.6 was determined. 
Unfortunately due to the dependence of mass-loss on the stel- 
lar luminosity and mass, a comparison of the known Galactic 
B-type supergiant mass-loss rates and those derived here will 
not depend solely on their different metallicities. 

Vink et al. (2001) provide a recipe to determine the mass- 
loss rate which includes a metallicity dependence of 0.85. 
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Fig. 10. To show the effect of incoherent electron scattering on the wings of AV362 and AV22. Shown are the best fits (solid line; 

) using the parameters presented in Table[3]&|4]and the same model with the inclusion of electron scattering (dashed line; - 

- -). Note the improvement in the fit to the wings of AV22 by including electron scattering and the little effect the scattering has 
on the fit to AV362. 



Three different formulae are used to describe the stellar mass- 
loss depending on its position relative to the so called 'bi- 
stability jumps', which vary due to the temperature and den- 
sity regime of the star. In Fig. [8] we compare our observed 
mass-loss rates with these theoretical predictions. Note that for 
the early B-type stars (BO - Bl) Vink's predictions, which are 
based on the recipe for stars above the first jump, underesti- 
mates M by a factor of 4, whilst for the mid B-type stars the 
recipe for between the two bistability jumps overestimates M 
by a factor of 10. This underestimation for the early B-type 
supergiants has also been observed for O-type supergiants and 
dwarfs by Bouret et al. (2003) and Evans et al. (2003). If we 
use the steeper exponent of 0.5 for the metallicity dependence 
in Vink's recipes we get good agreement with the mass-loss 
rates for the early B-type supergiants but the rates for the mid 
B-type supergiants are now over estimated by a factor of 16. 
Clearly, a better understanding of the metallicity and temper- 
ature dependence of mass-loss rates for B-type supergiants is 
required. 

An alternative perspective is provided by considering the 
wind momentum, D mom = (Mv^R 1 / 2 ). As described in the in- 
troduction, radiative line-driven wind theory relates the wind- 
momentum of a star to its luminosity. In recent years massive 
star analyses have endeavoured to calibrate this relationship as 
a function of spectral type and metallicity, with the plan to in- 
voke it as a method of deriving distances from OBA-type su- 
pergiants out to the Virgo & Fornax clusters (Puis et al. 1996; 
Kudritzki et al. 1999, 2000; Crowther et al. 2002; Bouret et al. 
2003; Bresolin et al. 2001; Urbaneja et al. 2003; Evans et al. 
2003). In Fig. [9] we compare the wind-momenta of the SMC 
B-type supergiants/giant studied here with the Galactic anal- 
ysis of Kudritzki et al. (1999). The solid lines in this figure 
represent the Galactic WLR relationships as derived from the 
recent line-blanketed O-type star analysis of Repolust, Puis, & 



Herrero (2003) and the unblanketed early and mid B-type su- 
pergiant analyses of Kudritzki et al. (1999, 2000). Surprisingly, 
the SMC stars from B0 - B2 agree well with the galactic re- 
lationship. However the two cool stars (B3 - B5) have wind- 
momenta which lie above the mid B relationship, this is more 
significant in the case of the B5 star, AV22. We must note that 
the mid B WLR of Kudritzki et al. (1999) was based on spec- 
tral types no later than B3, with the coolest star in their sam- 
ple, HD 14134, being at 18 kK in comparison to the 14 kK 
of AV362. From Fig.|9]we see that AV22 agrees well with the 
WLR calibrated by Kudritzki et al. (1999) for A-type stars. It is 
plausible that this cooler spectral type should follow the A-type 
star WLR or indeed a different calibration of the WLR for late 
B-type stars. Yet without more stars of similar or cooler spec- 
tral types (B5-B9), we can not make any definite conclusions. 
Moreover the inclusion of clumping into our analysis would 
decrease the mass-loss rate in such a way as to bring the wind- 
momenta of AV362 and AV22 into better agreement with the 
mid B WLR. 

From Fig.|3]it is apparent that the (far) wings of H a in these 
stars are quite strong and are unaccounted for by the FASTWIND 
models. In Fig. ^| we show the results of including incoherent 
electron scattering in the FASTWIND calculations. The electron 
scattering only increases the strength of the continuum around 
the H Q by 2% and 1% in AV362 and AV22, respectively, where 
as the observed spectrum shows a raise in the H Q wings of 7% 
and 3%. Whilst the situation has improved for AV22, the elec- 
tron scattering does little to improve the fit to H Q in AV362. 
Note that the fact that the electron-scattering can not account 
for the full strength of the electron scattering wings most prob- 
ably is not related to any presence or non-presence of wind- 
clumping, since these wings are formed in the upper photo- 
sphere (cf. Kudritzki et al. 1999). 
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Fig. 9. Wind momentum as a function of Luminosity. The solid 
lines are the wind-momentum-luminosity relationship, in de- 
creasing values of wind momentum, from galactic O-type su- 
pergiants without clumping and with line-blanketing (Repolust, 
Puis, & Herrero 2003) and from galactic early and mid B-type 

supergiants (Kudritzki & Puis 2000). The dashed line ( ) 

represents the unblanketed WLR for O-type supergiants deter- 
mined by Puis et al. (1996), whilst the dotted line (• • •) repre- 
sents that for the A-type supergiants determined by Kudritzki 
et al. (1999). Note the galactic BA-type WLR's were calibrated 
using an unblanketed version of FASTWIND. The plot shows 
the early and mid B-type stars of the galaxy (o, unfilled stars ; 
Kudritzki & Puis 2000) and the SMC (• , * ; this work). 



The agreement of our analysis with the galactic WLR for 
early B-type stars maybe somewhat misleading as the analy- 
sis of Kudritzki et al. (1999) was undertaken using an unblan- 
keted version of FASTWIND. Repolust, Puis, & Herrero (2003) 
and Puis et al. (1996) have carried out a line-blanketed and un- 
blanketed analysis, respectively, of a set of Galactic O stars. 
In Fig. |9] the regressions fits of both analyses are plotted and 
imply that the line-blanketed version of the O-type star WLR 
is steeper. Additionally for a given luminosity, higher wind- 
momenta are required than for the unblanketed results. Since a 
similar effect on the WLR is expected for the B-type stars, it is 
therefore possible that this is masking any predicted metallic- 
ity effect. In addition the uncertainty in the adopted distance to 
the SMC introduces a systematic error to the derived radii and 
mass-loss rates, this will also effect the wind-momentum lumi- 
nosity relationship. This changes the WLR through the wind- 
momentum and the luminosity, shifting the relationship along 
both the horizontal and vertical directions of Fig. [9] 

In general, the fitting uncertainties on the mass-loss rates 
are of the order 10 - 15 % , however the assumption of no 
clumping in our models introduces an additional uncertainty. 
Crowther et al. (2002) and Evans et al. (2003) noticed some in- 
dication of clumping in their O- and B-type supergiants from 
observations in the UV and far UV (see also Massa et al. 
2003, Hillier et al. 2003). Moreover Repolust, Puis & Herrero 



(2003) suggested the presence of clumping to explain certain 
inconsistencies between observed and theoretically predicted 
WLR's. Since the mass-loss rate varies with the filling factor as 
M/^/f , assuming a filling factor of 0. 1 reduces the mass-loss to 
one-third of its value before considering clumping. Although 
clumping is an important factor to be considered in the anal- 
ysis of supergiants with thick winds, the degree of clumping 
is still somewhat uncertain. Whilst for WR stars the clumping 
enhances the strength of the coherent electron scattered wings 
of the lines, this effect (which must not be confused with the 
incoherent electron scattering wings discussed above) is too 
small for detection in the OB-type supergiants (see Hillier et 
al. 2003). In addition by including clumping into the analyses 
of massive stars the mass-loss rate will be reduced and hence 
we would expect a shift in the WLR to lower wind-momenta 
(Repolust, Puis, & Herrero 2003). 

7. Conclusions 

From high-resolution UVES data we have derived the atmo- 
spheric and wind parameters of seven SMC B-type supergiants 
and one giant in the spectral range B0 - B5. In addition we have 
carried out an abundance analysis on each of these stars. This 
analysis was undertaken using the non-LTE line-blanketed at- 
mospheric code FASTWIND, which through spherical symme- 
try and extended atmospheres allows for the computation of the 
stellar wind (Santoloya-Rey et al. 1997; Herrero et al. 2002; 
Repolust, Puis, & Herrero 2003). 

The absolute abundances of the a-processed elements (O, 
Mg & Si) and carbon in these SMC stars are in good agreement 
with SMC H II regions and an unprocessed B-type star, AV304 
(Kurt et al. 1999; Rolleston et al. 2003). The exception is nitro- 
gen which in the giant, AV216, is enhanced by a factor of four 
while there is a mean overabundance of a factor of 13 present 
in the supergiants. Slightly lower nitrogen enhancements are 
observed in B-type giants in NGC330 and A-type supergiants 
(Lennon et al. 2003; Venn et al. 1999). 

A comparison was carried out with the most recent stellar 
evolution models of Maeder & Meynet (2001) which include 
rotationally induced mixing. However we recalibrated these 
models to an SMC initial nitrogen abundance of 6.55 dex as ob- 
served in AV304 instead of the one-fifth solar implemented by 
Maeder & Meynet. These models can for the most part repro- 
duce the observed nitrogen abundances if an initial rotational 
velocity of 300 kms -1 is invoked. Nevertheless, observed pro- 
jected rotational velocities are < 100 kms -1 and therefore we 
conclude that these models need to have more efficient mixing 
at lower initial rotational velocities. In addition we note that 
the evolutionary masses are generally higher than the derived 
spectroscopic masses by at least a factor of 1.2. Some of this 
discrepancy may be accounted for in the errors of the adopted 
distance modulus and possibly the gravity, leaving a residual 
difference in masses of less than 0. 1 dex where M evo i > M spec . 

The theoretical mass-loss rate predictions of Vink et al. 
(2001) which implement a metallicity dependence of Z° 64 , un- 
derestimate M for the early B-type stars and overestimate it 
for the mid B-types. By using a steeper metallicity dependence 
(Z 5 ), the early B-type stars are in agreement but the situation 
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is then worse for the later spectral types. Clearly the metallicity 
dependence of mass-loss rates for B-type stars is still uncertain 
and a comparison of equivalent sets of stars in different metal- 
licity environments is needed to understand this dependence. 

Finally, we have compared our results to the wind- 
momentum relationships of Galactic B-type stars by Kudritzki 
et al. (1999, 2000). The wind momenta of our stars are in good 
agreement with the galactic B-type supergiants. This is a puz- 
zling result given the theoretical predictions of a metallicity 
dependence on the wind-momenta and the observational evi- 
dence of this effect for O-type supergiants in the Magellanic 
Clouds (Puis et al. 1996; Bouret et al. 2003). However the ob- 
served Galactic WLR was derived using an unblanketed anal- 
ysis which may mask any predicted metallicity dependence 
when compared to our results. 

Work is currently underway to analyse EMMI spectra of 
additional SMC early B-type supergiants in the spectral range 
B0 - B2.5. This larger sample of stars will hopefully provide 
more insight into the behaviour of the wind-momenta and stel- 
lar evolution of these early B-type supergiants at low metallic - 
ities. 
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